program test
use general
implicit none
integer,parameter :: N=16
double precision :: pi
double precision :: xp,length,dx,dy,yp
integer :: xgrid,xgridn,ip,icounter,ygrid,ygridn
double precision :: ddx,ddy
double precision, parameter :: n1=16,n2=16
double precision, dimension(n1+2,n2) :: ux,uy,omega
double precision, dimension(3) :: A,B,C
A=0.;B=0.;C=0.
C=cross(A,B)
write(*,*)'A',A
write(*,*)'B',B
write(*,*)'C',C
endprogram test
!!-----------------------------------------------------
